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Abstract 

BDDC method is the most advanced method from the Balancing family of iterative 
substructuring methods for the solution of large systems of linear algebraic equations arising 
from discretization of elliptic boundary value problems. In the case of many substructures, 
solving the coarse problem exactly becomes a bottleneck. Since the coarse problem in BDDC 
has the same structure as the original problem, it is straightforward to apply the BDDC method 
recursively to solve the coarse problem only approximately. In this paper, we formulate a new 
family of abstract Multispace BDDC methods and give condition number bounds from the 
abstract additive Schwarz preconditioning theory. The Multilevel BDDC is then treated as 
a special case of the Multispace BDDC and abstract multilevel condition number bounds are 
given. The abstract bounds yield polylogarithmic condition number bounds for an arbitrary 
fixed number of levels and scalar elliptic problems discretized by finite elements in two and 
three spatial dimensions. Numerical experiments confirm the theory. 

AMS Subject Classification: 65N55, 65M55, 65Y05 

Key words: Iterative substructuring, additive Schwarz method, balancing domain 
decomposition, BDD, BDDC, Multispace BDDC, Multilevel BDDC 



1 Introduction 

The BDDC (Balancing Domain Decomposition by Constraints) method by Dohrmann [4] is the 
most advanced method from the BDD family introduced by Mandel [12] . It is a Neumann-Neumann 
iterative substructuring method of Schwarz type [5] that iterates on the system of primal variables 
reduced to the interfaces between the substructures. The BDDC method is closely related to 
the FETI-DP method (Finite Element Tearing and Interconnecting - Dual, Primal) by Farhat et 
al. [6l [7]. FETI-DP is a dual method that iterates on a system for Lagrange multipliers that 
enforce continuity on the interfaces, with some "coarse" variables treated as primal, and it is 
a further development of the FETI method by Farhat and Roux [8]. Polylogarithmic condition 
number estimates for BDDC were obtained in \13\ [T4"] and a proof that the eigenvalues of BDDC 
and FETI-DP are actually the same except for eigenvalue equal to one was given in Mandel et 
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al. [Hj. Simpler proofs of the equality of eigenvalues were obtained by Li and Widlund [10], and 
also by Brenner and Sung [3], who also gave an example when BDDC has an eigenvalue equal 
to one but FETI-DP does not. In the case of many substructures, solving the coarse problem 
exactly becomes a bottleneck. However, since the coarse problem in BDDC has the same form as 
the original problem, the BDDC method can be applied recursively to solve the coarse problem 
only approximately. This leads to a multilevel form of BDDC in a straightforward manner, see 
Dohrmann [3]. Polylogarithmic condition number bounds for three-level BDDC (BDDC with two 
coarse levels) were proved in two and three spatial dimensions by Tu (20l [19] . 

In this paper, we present a new abstract Multispace BDDC method. The method extends a 
simple variational setting of BDDC from Mandel and Sousedfk |15j . which could be understood as 
an abstract version of BDDC by partial subassembly in Li and Widlund [llj. However, we do not 
adopt their change of variables, which does not seem to be suitable in our abstract setting. We 
provide a condition number estimate for the abstract Multispace BDDC method, which generalizes 
the estimate for a single space from |15] , The proof is based on the abstract additive Schwarz 
theory by Dryja and Widlund [5]. Many BDDC formulations (with an explicit treatment of 
substructure interiors, after reduction to substructure interfaces, with two levels, and multilevel) 
are then viewed as abstract Multispace BDDC with a suitable choice of spaces and operators, and 
abstract condition number estimates for those BDDC methods follow. This result in turn gives 
a polylogarithmic condition number bound for Multilevel BDDC applied to a second-order scalar 
elliptic model problems, with an arbitrary number of levels. A brief presentation of the main results 
of the paper, without proofs and with a simplified formulation of the Multispace BDDC estimate, 
is contained in the conference paper [16]. 

The paper is organized as follows. In Sec. [2] we introduce the abstract problem setting. In Sec. [3] 
we formulate an abstract Multispace BDDC as an additive Schwarz preconditioner. In Sec. H] we 
introduce the settings of a model problem using finite element discretization. In Sec. [5] we recall 
the algorithm of the original (two- level) BDDC method and formulate it as a Multispace BDDC. 
In Sec. [6] we generalize the algorithm to obtain Multilevel BDDC and we also give an abstract 
condition number bound. In Sec. [7] we derive the condition number bound for the model problem. 
Finally, in Sec. [HJ we report on numerical results. 

2 Abstract Problem Setting 

We wish to solve an abstract linear problem 

u G X : a(u,v) = (f x ,v) , Vi> G X, (1) 

where X is a finite dimensional linear space, a (•, •) is a symmetric positive definite bilinear form 
defined on X, fx & X' is the right-hand side with X' denoting the dual space of X, and (•, •) is the 
duality pairing. The form a (•, •) is also called the energy inner product, the value of the quadratic 

1 /2 

form a{u,u) is called the energy of u, and the norm ||u|| a = a{u,u) 1 is called the energy norm. 
The operator Ax '■ X i— > X' associated with a is defined by 

a(u,v) = (Axu,v) , Vu, v G X. 

A preconditioner is a mapping B : X' — > X and we will look for preconditioners such that 
(r, Br) is also symmetric and positive definite on X' . It is well known that then BAx '■ X — > X has 
only real positive eigenvalues, and convergence of the preconditioned conjugate gradients method 
is bounded using the condition number 

_ ^max(BAx) 
^min(BAx) ' 
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which we wish to bound above. 

All abstract spaces in this paper are finite dimensional linear spaces and we make no distinction 
between a linear operator and its matrix. 



3 Abstract Multispace BDDC 

To introduce abstract Multispace BDDC preconditioner, suppose that the bilinear form a is defined 
and symmetric positive semidefinite on some larger space W D X. The preconditioner is derived 
from the abstract additive Schwarz theory, however we decompose some space between X and 
W rather than X as it would be done in the additive Schwarz method: In the design of the 
preconditioner, we choose spaces V k , k = 1, . . . , M, such that 

M 



Xc^FfeCW. (2) 



k=\ 

Assumption 1 The form a(-, •) is positive definite on each V k separately. 

Algorithm 2 (Abstract Multispace BDDC) Given spaces Vk and linear operators Qk, k 
1, . . . , M such that a (•, •) is positive definite on each space Vk, and 

M 



xcJ2 v k, Qk:V k ^x, 



fe=i 

define the preconditioner B : r E X' \ — > u € X by 

M 



B :r^u = ^2 QkVk, v k G V k : a (v k , z k ) = (r, Q k z k ) , \/z k G V k . (3) 



k=l 



We formulate the condition number bound first in the full strength allowed by the proof. The 
bound used in the rest of this paper will be a corollary. 

Theorem 3 Define for all k = 1, . . . , M the spaces by 

V k M = { v k G V k : Vzfe G V k : Q k v k = Qk^k =^ < Ikfclla} • 

If there exist constants Co, u, and a symmetric matrix £ = (eij)f^ =1 , such that 

M M 

VueX lv k eV k , k = l,...,M:u = ^2Q k v k , ^2\\v k \\l<C \\u\\l (4) 

k=l k=l 

\fk = l,...,M Vv k eV k M : \\Q k v k \\ 2 a <u\\v k \\l (5) 
Vz k €Q k V k , k = l,...,M :a(zi,Zj) <eij\\zi\\ a \\zj\\ a , (6) 



then the preconditioner from Algorithmic satisfies 



< C iop(£). 
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Proof. We interpret the Multispace BDDC preconditioner as an abstract additive Schwarz method. 
An abstract additive Schwarz method is specified by a decomposition of the space X into subspaces, 

X = X 1 + ... + X M , (7) 

and by symmetric positive definite bilinear forms bi on Xj. The preconditioner is a linear operator 

B : X' -> X, B :r^u, 

defined by solving the following variational problems on the subspaces and adding the results, 

M 

B : r *—>■ u = ^ u k , u k G X k : b k (u k ,y k ) = (r, y k ) , \/y k G X k . (8) 
k=l 

Dryja and Widlund [5] proved that if there exist constants Co, u, and a symmetric matrix 
£ = (eij)fj = i, such that 

M M 

\/u G X 3u k G X k , k = 1, ... ,M : u = ^Uk, ^ IKI|& fc < Cd \\u\\ 2 a (9) 

fc=l fc=l 

Vk = 1, . . . ,M Vn fc G X k : \\uk\\l < lo \\v,k\\l h (10) 
Vufc G X k , k = l,...,M : a(ui,Uj) <eij\\ui\\ a \\uj\\ a (11) 

then 

Amax(-B^x) , n ( c\ 

where p is the spectral radius. 

Now the idea of the proof is essentially to map the assumptions of the abstract additive Schwarz 
estimate from the decomposition ([7]) of the space X to the decomposition (|2|). Define the spaces 

X k = Q k V k . 

We will show that the preconditioner ([3]) satisfies ([S]), where b k is defined by 

h(u k ,y k ) = a(G k x,G k z) , x,z£X, u k = Q k G k x, y k = Q k G k z. (12) 
with the operators G k : X — > V k defined by 

1 M 

G k :u^ v k , -a (v k , v k ) -> min, s.t. v k G V^, u = ^ QkVk, (13) 

k=l 

First, from the definition of operators G k , spaces X k , and because a is positive definite on V k by 
Assumption (JTJ), it follows that G k x and G k z in (fT2|) exist and are unique, so b k is defined correctly. 
To prove ([5D, let v k be as in ([3]) and note that v k is the solution of 

-a (v k , v h ) - (r, Q k v k ) -» min, w fc G 

Consequently, the preconditioner ([3]) is an abstract additive Schwarz method and we only need to 
verify the inequalities ([9^-(fTT]). To prove ©, let u G X. Then, with v k from the assumption Q 
and with u k = Q k G k v k as in (|12p . it follows that 

Af Af Af 

u = j^u fc , \\ u k\\b k = ii u *iia - c ° h2- 

fe=l fc=l k=l 
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Next, let Uk G X}~. From the definitions of Xj- and V^, it follows that there exist unique Vk € V^ 
such that Uk = QkVk- Using the assumption ([5]) and the definition of bk in (fT2|) . we get 

||«*|la = WQkVkWl < V \\vk\\l = UJ \\u k \\l k , 



which gives (jlOp . Finally, © is the same as (jlip . ■ 

The next Corollary was given without proof in |X6|, Lemma 1]. This is the special case of 
Theorem [3] that will be actually used. In the case when M = 1, this result was proved in [15j . 

Corollary 4 Assume that the subspaces Vk are energy orthogonal, the operators Qk are projections, 
a (•, •) is positive definite on each space Vk, and 



Vu£l: 



M 



u 



^Vk, v k G V k 



k=l 



M 



(14) 



k=l 



Then the abstract Multispace BDDC preconditioner from Algorithmic satisfies 



Amax(-B^x) 
^min{BAx) 



< CO 



max sup 

k v k ev k 



ii ii- 
\\m a 



(15) 



Proof. We only need to verify the assumptions of Theorem [3j Let u E X and choose Vk as 
the energy orthogonal projections of u on Vk- First, since the spaces Vk are energy orthogonal, 
u = Yl v ki Qk are projections, and from (fT4"l) u = J2Qk v k, we get that = X^lbfclla which 

proves dH) with Cq = 1. Next, the assumption © becomes the definition of u in (|15p . Finally, ([6]) 
with 6 = 1 follows from the orthogonality of subspaces Vk- ■ 

Remark 5 The assumption {1$ can be written as 

M 



Y,QkPk 



k=l 



X 



where Pk is the a-orthogonal projection from ®jLiVj onto Vk- Hence, the property {1$ is a type 
of decomposition of unity. 

In the case when M = 1, HIX) means that the projection Q\ is onto X . 



4 Finite Element Problem Setting 

Let n be a bounded domain in M. d , d = 2 or 3, decomposed into N nonoverlapping subdomains 
n s , s = 1, .--,N, which form a conforming triangulation of the domain $7. Subdomains will be also 
called substructures. Each substructure is a union of Lagrangian PI or Ql finite elements with 
characteristic mesh size h, and the nodes of the finite elements between substructures coincide. The 
nodes contained in the intersection of at least two substructures are called boundary nodes. The 
union of all boundary nodes is called the interface T. The interface T is a union of three different 
types of open sets: faces, edges, and vertices. The substructure vertices will be also called corners. 
For the case of regular substructures, such as cubes or tetrahedrons, we can use standard geometric 
definition of faces, edges, and vertices; cf., e.g., [9] for a more general definition. 

In this paper, we find it more convenient to use the notation of abstract linear spaces and linear 
operators between them instead of the space 1" and matrices. The results can be easily converted 
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to the matrix language by choosing a finite element basis. The space of the finite element functions 
on f2 will be denoted as U. Let W s be the space of finite element functions on substructure Q s , 
such that all of their degrees of freedom on d£l s H dO, are zero. Let 

W = W 1 x ••• x W N , 

and consider a bilinear form arising from the second-order scalar elliptic problem as 

a(u,v) = y / VuVvdx, u, v £ W. (16) 

Now U C W is the subspace of all functions from W that are continuous across the substructure 
interfaces. We are interested in the solution of the problem (pQ) with X = U, 

u<EU :a(u,v) = (f,v) , W G U, (17) 

where the bilinear form a is associated on the space U with the system operator A, defined by 

A : U !->• U', a(u, v) = (Au, v) for all u,v e U, (18) 

and / € U' is the right-hand side. Hence, (fTTl) is equivalent to 

An = /. (19) 

Define Ui C U as the subspace of functions that are zero on the interface T, i.e., the "interior" 
functions. Denote by P the energy orthogonal projection from W onto Uj, 

P : w £ W i — ► vi £ Ui : a (vi, zf) = a (w, Zi) , Vzj € f/j. 

Functions from (/ — P) W, i.e., from the nullspace of P, are called discrete harmonic; these functions 
are a-orthogonal to Uj and energy minimal with respect to increments in Uj. Next, let W be the 
space of all discrete harmonic functions that are continuous across substructure boundaries, that is 

W = (I - P) U. (20) 

In particular, 

U = U T ®W, Hi ± a W. (21) 

A common approach in substructuring is to reduce the problem to the interface. The problem 
(fTT|) is equivalent to two independent problems on the energy orthogonal subspaces Ui and W, and 
the solution u satisfies u = uj + u, where 



u^Ur. a( Ul , Vl ) = {f, VI ) , Vu/ G U h (22) 
u G W : a(u,v) = (f,v) , VveW. (23) 

The solution of the interior problem (I22jl decomposes into independent problems, one per each 
substructure. The reduced problem ([23]) is then solved by preconditioned conjugate gradients. The 
reduced problem (|23j) is usually written equivalently as 

u G W : s(u, v) = (g, v ) , \/v G W, 
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where s is the form a restricted on the subspace W, and g is the reduced right hand side, i.e., the 
functional / restricted to the space W. The reduced right-hand side g is usually written as 

(g,v) = (f,v)-a( Ul ,v), VveW, (24) 

because a(ui, v) = by (|2ip . In the implementation, the process of passing to the reduced problem 
becomes the elimination of the internal degrees of freedom of the substructures, also known as 
static condensation. The matrix of the reduced bilinear form s in the basis defined by interface 
degrees of freedom becomes the Schur complement, and (|24|) becomes the reduced right-hand side. 
For details on the matrix formulation, see, e.g., [TTJ Sec. 4.6] or [181 Sec. 4.3]. 

The BDDC method is a two-level preconditioner characterized by the selection of certain coarse 
degrees of freedom, such as values at the corners and averages over edges or faces of substructures. 
Define W C W as the subspace of all functions such that the values of any coarse degrees of freedom 
have a common value for all relevant substructures and vanish on <9f2, and Wa C W as the subspace 
of all function such that their coarse degrees of freedom vanish. Next, define Wn as the subspace 
of all functions such that their coarse degrees of freedom between adjacent substructures coincide, 
and such that their energy is minimal. Clearly, functions in Wn are uniquely determined by the 
values of their coarse degrees of freedom, and 

W A ^-a W n , and W = W A W n . (25) 

We assume that 

a is positive definite on W. (26) 

That is the case when a is positive definite on the space U, where the problem ([I]) is posed, and 
there are sufficiently many coarse degrees of freedom. We further assume that the coarse degrees 
of freedom are zero on all functions from J7j, that is, 

Ul C W A . (27) 

In other words, the coarse degrees of freedom depend on the values on substructure boundaries 
only. From (|25ft and (|27p . it follows that the functions in Wu are discrete harmonic, that is, 

W n = (I-P)W n . (28) 

Next, let E be a projection from W onto U, defined by taking some weighted average on 
substructure interfaces. That is, we assume that 

E : W -»• U, EU = U, E 2 = E. (29) 

Since a projection is the identity on its range, it follows that E does not change the interior degrees 
of freedom, 

EUi = U T , (30) 

since Ui C U. Finally, we show that the operator (/ — P) E is a projection. From (|30p it follows 
that E does not change interior degrees of freedom, so EP = P. Then, using the fact that I — P 
and E are projections, we get 

[(I - P) E} 2 = (I - P) E (I — P) E 

= {I-P){E-P)E (31) 
= (I - P)(I - P)E = (I - P) E. 
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Remark 6 In \14\ UU!f , the whole analysis was done in spaces of discrete harmonic functions after 
eliminating Uj, and the space W was the solution space. In particular, W consisted of discrete 
harmonic functions only, while the same space here would be (I — P)W. The decomposition of this 
space used in \T1\ [151/ would be in our context written as 

(/ - P)W = (/ - P)W A Wn, (/ - P)W A ± a W n . (32) 

In the next section, the space X will be either U or W. 

5 Two-level BDDC as Multispace BDDC 

We show several different ways the original, two-level, BDDC algorithm can be interpreted as 
multispace BDDC. We consider first BDDC applied to the reduced problem (f23]) . that is, (P) with 
X = W. This was the formulation considered in [14J. Define the space of discrete harmonic 
functions with coarse degrees of freedom continuous across the interface 

W V = {I- P) W. 

Because we work in the space of discrete harmonic functions and the output of the averaging 
operator E is not discrete harmonic, denote 

E V = {I- P) E. (33) 

In an implementation, discrete harmonic functions are represented by the values of their degrees of 
freedom on substructure interfaces, cf., e.g. [18] ; hence, the definition ([33]) serves formal purposes 
only, so that everything can be written in terms of discrete harmonic functions without passing to 
the matrix formulation. 

Algorithm 7 (|15j, BDDC on the reduced problem) Define the preconditioner r 6 W' i — ► 
u £ W by 

u = Epwr, wr £ Wr '■ a (wr, zr) = (r, E^zr) , Vzp € Wr- (34) 

Proposition 8 ( [15J ) The BDDC preconditioner on the reduced problem in Algorithm^ is the 
abstract Multispace BDDC from Algorithm^ with M = 1 and the space and operator given by 

X = W, Vx = W r , Qi = E T . (35) 

Also, the assumptions of Corollary are satisfied. 

Proof. We only need to note that the bilinear form a(-, •) is positive definite on Wr C W by (j26|) . 
and the operator Er defined by ([33]) is a projection by (f3~T[) . The projection Er is onto W because 
E is onto U by ([29]) . and I — P maps U onto W by the definition of W in ([201) . ■ 

Using the decomposition ([32]) . we can split the solution in the space Wr into the independent 
solution of two subproblems: mutually independent problems on substructures as the solution in 
the space Wta = (I — P)W&, and a solution of global coarse problem in the space Wn- The space 
Wr has a decomposition 

W t = Wta®Wu, and W rA ± a Wn, (36) 
the same as the decomposition (132j) . and Algorithm [7] can be rewritten as follows. 
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Algorithm 9 (|14j, BDDC on the reduced problem) Define the preconditioner r G W i — > 
u & W by u = Er (wrA + ^n); where 

wrA G WrA ■ a (wrA, zfa) = {r, -Et^ta) , VzrA G WrA, (37) 
i«n G Wu ■ a (w n , zrn) = (r, E r z r n) , Vz r n G W n - (38) 

Proposition 10 T/ie BDDC preconditioner on the reduced problem in Algorithm\Mis the abstract 
Multispace BDDC from Algorithm^ with M = 2 and the spaces and operators given by 

X = W, Vi = W T a, V 2 = Wn, Qi = Q2 = E v . (39) 
Also, the assumptions of Corollary^ are satisfied. 

Proof. Let r G W' . Define the vectors Vj, i = 1,2 in Multispace BDDC by ([3]) with Vi and Qi 
given by ([39]) . Let u, u>rA> wn be the quantities in Algorithm [9l defined by ([37]) - ([38]) . Using the 
decomposition ([36]) . any tor G Wr can be written uniquely as wr = wrA+wu for some wrA and 
wu corresponding to ([3]) as v\ = wrA and V2 = wn, and u = E? (wrA + wn)- 

To verify the assumptions of Corollary U note that the decomposition (|36p is a— orthogonal, 
o(-, •) is positive definite on both Wta and Wu as subspaces of Wr by ([26]) . and Ep is a projection 
by ([HI]). ■ 

Next, we present a BDDC formulation on the space U with explicit treatment of interior 
functions in the space Uj as in [HE], i- e -j i n the way the BDDC algorithm was originally formulated. 

Algorithm 11 ([4, 13J, original BDDC) Define the preconditioner r G V 1 — > u £ U as 

follows. Compute the interior pre- correction 

uj G £/> : a z z ) = (r, z z ) , Vzj G C/j. (40) 

iSet up t/ie updated residual 

rB^U', (tb,v) = (r,v) — a(ui,v) , V« £ 17. (41) 
Compute the substructure correction 

ua = Ewa, wa G Wa : a (wa, za) = (?*B) ^a) , Vza G Wa- (42) 
Compute the coarse correction 

uu = Ewu, wu G Wu ■ a (w n , zn) = (rs, #z n ) , Vz n G Wu- (43) 
^Idd i/ie corrections 

ub = ua + ^n- 

Compute the interior post-correction 

vi eUi : a(vi,zi) = a(u B ,zi) , Vzj G £//. (44) 
Apply the combined corrections 

u = ub — vj + itj. (45) 
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The interior corrections (|40p and (|44p decompose into independent Dirichlet problems, one 
for each substructure. The substructure correction (j42]) decomposes into independent constrained 
Neumann problems, one for each substructure. Thus, the evaluation of the preconditioner requires 
three problems to be solved in each substructure, plus solution of the coarse problem (1431) . In 
addition, the substructure corrections can be solved in parallel with the coarse problem. 

Remark 12 As it is well known $4%, the first interior correction ( L/0| ) can be omitted in the 
implementation by starting the iterations from an initial solution such that the residual in the 
interior of the substructures is zero, 

a(u,zi) - (fx,zi) = 0, Vzj £ C/j, 

i.e., such that the error is discrete harmonic. Then the output of the preconditioner is discrete 
harmonic and thus the errors in all the CG iterations (which are linear combinations of the original 
error and outputs from the preconditioner) are also discrete harmonic by induction. 

The following proposition will be the starting point for the multilevel case. 

Proposition 13 The original BDDC preconditioner in Algorithm [77] is the abstract Multispace 
BDDC from Algorithm^ with M = 3 and the spaces and operators given by 

X = U, V 1 = U I , V 2 = (I-P)W A , V 3 = W U , (46) 
Qi = 1, Q2 = Qs = (I- P) E, (47) 

and the assumptions of Corollary [|] are satisfied. 

Proof. Let r € U'. Define the vectors Vi, i = 1,2,3, in Multispace BDDC by ([3]) with the spaces 
Vi given by ([4"6]) and with the operators Qi given by (|4"T|). Let ui, rs, wa, wjj, ub, vi, and u be 
the quantities in Algorithm [UJ defined by (l4*0l) - (l4"5j) . 

First, with V\ = Z7j, the definition of v\ in ([3j) with k = 1 is identical to the definition of uj in 
(f40l) . so ui = v\. 

Next, consider wa € Wa defined in (fl^l) . We show that wa satisfies ([3]) with k = 2, i.e., 
V2 = wa- So, let za G Wa be arbitrary. From (|4"2"|) and (JUJ), 

a {wa,za) = {tb, Eza) = (r, Eza) - a (ui, Eza) ■ (48) 

Now from the definition of ui by (|40p and the fact that PEza G Ui, we get 

(r, P£z A ) - a ( UI , PEza) = 0, (49) 

and subtracting ()49p from (|48p gives 

a (u; A , z a ) = (r, (/ - P) Eza) - a (m, (I - P) Ez A ) 
= (r, (/ - P) Eza) , 

because a (uj, (I — P) Eza) = by orthogonality. To verify ([3]), it is enough to show that Pwa = 0; 
then wa S (J — P)Wa = Vi- Since P is an a-orthogonal projection, it holds that 

a (Pw A , Pwa) = a (wa, Pwa) = (r B , EPw A ) = 0, (50) 
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where we have used EUi C Uj following the assumption (|30p and the equality 

(tb, zi) = (r, zi) - a (ui, Zj) = 

for any zi € Uj, which follows from (|4ip and (|40p . Since a is positive definite on W D 17/ by 
assumption ([26]) . it follows from ([50]) that Pj«a = 0. 

In exactly the same way, from (li"3j) - (f4"5j) . we get that if »n G Wn is defined by ([431) . then 
^3 = "^n satisfies ([3]) with fc = 3. (The proof that Pwn = can be simplified but there is nothing 
wrong with proceeding exactly as for wa.) 

Finally, from ([44"]) , vi = P (Pwa + Ewu ) , so 

u = ui + [ub - vi) 
= ui + (I — P) Pw A + (I-P) Ew n 

= QlV! + Q 2 V 2 + Q3V3- 

It remains to verify the assumptions of Corollary SJ 

First, the spaces Wjj and Wa are a-orthogonal by ([25]) and, from ([27]) . 

thus (I — P) Wa J-a Wq. Clearly, (J — P) Wa J-a 17/. Since Wn consists of discrete harmonic 
functions from ([28]) . so Wn -L a ?7j, it follows that the spaces Vi, i = 1,2,3, given by ([4"6]) . are 
a-orthogonal. 

Next, (I — P) E is by (f3Tj) a projection, and so are the operators from ([4"T]) . 
It remains to prove the decomposition of unity (|14p . Let 

u^uz + WA + wnCt/, «/ G P/, w A G (I-P)Wa, w n G W n , (51) 

and let 

u = u/ + (I — P) £w A + (I - P) Pw n . 

From ([5T]) . wa + £ P since u' £ U and u/ € {7/ C P. Then P (wa + »n) = + by ([29]) . 

so 

u = uj + (I - P) E (w A + w n ) 
= «j + (I - P) (w A + w n ) 
= ui + wa + wn = u', 

because both wa and wn are discrete harmonic. ■ 

The next Theorem shows an equivalence of the three Algorithms introduced above. 

Theorem 14 The eigenvalues of the preconditioned operators from Algorithm^ and Algorithm^ 
are exactly the same. They are also the same as the eigenvalues from Alaorithm \lll except possibly 
for multiplicity of eigenvalue equal to one. 

Proof. From the decomposition (|36p . we can write any w G Wp uniquely as w = wa + wn for some 
wa G Wta and wn G Wn, so the preconditioned operators from Algorithms and are spectrally 
equivalent and we need only to show their spectral equivalence to the preconditioned operator from 
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Algorithm II li First, we note that the operator A : U 
form as 



V defined by (|18p . and given in the block 



with blocks 



A 

A n : Uz 
An : W 



A n A IT 
An Arr 



U'n 
U'j, 



Air ■ U x 
Arr ■ W 



W, 

w' 



is block diagonal and An = Air = for any u G U, written as u = ui + w, because Ui _L a W. 
Next, we note that the block ^4rr : W W is the Schur complement operator corresponding to 
the form s. Finally, since the block An is used only in the preprocessing step, the preconditioned 
operator from Algorithms [7] and [9] is simply MrrArr : r G W — > u G W. 

Let us now turn to Algorithm 1111 Let the residual r G U be written as r = rj + rr, where 
rj G U'j and rr € W. Taking rr = 0, we get r = rj, and it follows that r# = ub = vj = 0, 
so u = ui. On the other hand, taking r = rr gives ui = 0, r# = rr, vi = Pub and finally 
u = (I — P) E(wa + wjj), so u G W. This shows that the off-diagonal blocks of the preconditioner 
M are zero, and therefore it is block diagonal 

Mn 
M rr 



M 



Next, let us take u = uj, and consider rr = 0. The algorithm returns tb = ub 
u = uj. This means that MnAnui = uj, so Mn = Ajj 1 . The operator A : U 



preconditioned operator MA : r £ U' 

An 
^rr 



vi = 0, and finally 
U', and the block 

u G U from Algorithm [TT] can be written, respectively, as 



MA 



I 
MrrArr 



where the right lower block Mrr^rr : r G W — > u G W is exactly the same as the preconditioned 
operator from Algorithms [7] and [9j ■ 

The BDDC condition number estimate is well known from [13J. Following Theorem 1141 and 
Corollary U we only need to estimate — P) Ew\\ on W. 

Theorem 15 (|13j) The condition number of the original BDDC algorithm satisfies k <oo, where 



uJ 



max ■ 



sup ■ 



(I-P)Ew\\l 



1} <C 1 + log 



w 



H 

T 



(52) 



Remark 16 In J13\j . the theorem was formulated by taking the supremum over the space of discrete 
harmonic functions (I — P)W. However, the supremum remains the same by taking the larger space 
WD(I- P)W, since 

\\(I-P)Ew\\ 2 a < \\(I-P)E(I-P)w\\l 



\\W\\-a W-P)v\fa 

from E (I — P) = E, which follows from l3ty). and from \\w\\ a > \ 
the a-orthogonality of the projection P. 



(I — P) w\\ , which follows from 



Before proceeding into the Multilevel BDDC section, let us write concisely the spaces and 
operators involved in the two-level preconditioner as 

Ui^U^W A ®W n = W CW. 
We are now ready to extend this decomposition into the multilevel case. 
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6 Multilevel BDDC and an Abstract Bound 



In this section, we generalize the two-level BDDC preconditioner to multiple levels, using the 
abstract Multispace BDDC framework from Algorithm [2j The substructuring components from 
Section [5] will be denoted by an additional subscript i, as Of , s = 1, . . . N±, etc., and called level 1. 
The level 1 coarse problem ()43|) will be called the level 2 problem. It has the same finite element 
structure as the original problem ([T|) on level 1, so we put Ui = Wni- Level 1 substructures are 
level 2 elements and level 1 coarse degrees of freedom are level 2 degrees of freedom. Repeating this 
process recursively, level i — 1 substructures become level i elements, and the level i substructures 
are agglomerates of level i elements. Level i substructures are denoted by Qf, s = 1, . . . ,Ni, and 
they are assumed to form a conforming triangulation with a characteristic substructure size Hi. 
For convenience, we denote by S7g the original finite elements and put Hq = h. The interface Tj on 
level i is defined as the union of all level % boundary nodes, i.e., nodes shared by at least two level i 
substructures, and we note that Tj C IV-i. Level i — 1 coarse degrees of freedom become level i 
degrees of freedom. The shape functions on level i are determined by minimization of energy with 
respect to level i — \ shape functions, subject to the value of exactly one level i degree of freedom 
being one and others level i degrees of freedom being zero. The minimization is done on each level 
i element (level i — 1 substructure) separately, so the values of level i — 1 degrees of freedom are 
in general discontinuous between level i — 1 substructures, and only the values of level i degrees of 
freedom between neighboring level i elements coincide. 

The development of the spaces on level i now parallels the finite element setting in Section [H 
Denote f7, = Wm-i- Let Wf be the space of functions on the substructure Qf, such that all of 
their degrees of freedom on dQf n d£l are zero, and let 

Wi = W/ x ••• x Wf\ 

Then Ui C Wj is the subspace of all functions from W that are continuous across the interfaces Tj. 
Define Un C Ui as the subspace of functions that are zero on Tj, i.e., the functions "interior" to 
the level i substructures. Denote by Pi the energy orthogonal projection from Wj onto Un, 

Pi : Wi G Wi 1 — ► vn eUu:a (vu, z R ) = a (wi, Zu) , Vzn € Un. 

Functions from (/ — Pj) Wi, i.e., from the nullspace of Pi, are called discrete harmonic on level 
i) these functions are a-orthogonal to Un and energy minimal with respect to increments in Uu- 
Denote by W% C Ui the subspace of discrete harmonic functions on level i, that is 

Wi = (I- Pi) Ui. (53) 

In particular, Un _L a Wj. Define Wj C Wj as the subspace of all functions such that the values of 
any coarse degrees of freedom on level i have a common value for all relevant level i substructures 
and vanish on <9f2f n d£l, and Wai C Wj as the subspace of all functions such that their level i 
coarse degrees of freedom vanish. Define Wm as the subspace of all functions such that their level 
i coarse degrees of freedom between adjacent substructures coincide, and such that their energy is 
minimal. Clearly, functions in Wm are uniquely determined by the values of their level i coarse 
degrees of freedom, and 

W Al ± a W m , Wi = W Al Wm- (54) 
We assume that the level i coarse degrees of freedom are zero on all functions from Un, that is, 

U u c W Ai . (55) 
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Figure 1: Space decompositions, embeddings and projections in the Multilevel BDDC 

In other words, level i coarse degrees of freedom depend on the values on level i substructure 
boundaries only. From (|54p and (|55p . it follows that the functions in Wm are discrete harmonic on 
level i, that is 

Wm = {I-Pi)Wm- (56) 
Let E be a projection from Wi onto U{, defined by taking some weighted average on Tj 

Ei-.Wi^ U h EiUji = Uu, El = Ei. 

Since projection is the identity on its range, Ei does not change the level i interior degrees of 
freedom, in particular 

EiUn = Uu. (57) 

Finally, we introduce an interpolation Jj : Ui — > Ui from level i degrees of freedom to functions 
in some classical finite element space Ui with the same degrees of freedom as Ui. The space Ui will 
be used for comparison purposes, to invoke known inequalities for finite elements. A more detailed 
description of the properties of Ii and the spaces Ui is postponed to the next section. 

The hierarchy of spaces and operators is shown concisely in Figure HJ The Multilevel BDDC 
method is defined recursively [U [16] by solving the coarse problem on level i only approximately, by 
one application of the preconditioner on level i — 1. Eventually, at level, L — l, the coarse problem, 
which is the level L problem, is solved exactly. We need a more formal description of the method 
here, which is provided by the following algorithm. 

Algorithm 17 (Multilevel BDDC) Define the preconditioner n 6 U[ 1 — > u\ G U\ as follows: 
for i = l,...,L-l, 

Compute interior pre- correction on level i, 

uu eUn:a (un, zn) = (n, z R ) , \/z R e Uu. (58) 
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Get updated residual on level i, 

rBi G Ui, (r B i, Vi) = (n, Vi) - a (un, Vi) , Vt?* G U{. (59) 

Find i/ie substructure correction on level i: 

WAi G Wai : a (tUAi, ZAi) = (rBi, £^Ai) , VzAi G W^Ai- (60) 
Formulate the coarse problem on level i, 

wm G Wm ■ a (wm, z m ) = (rBi,EiZm) , Vz m G Wrii, (61) 

If i = L — 1, solve the coarse problem directly and set ul = wni-i, 
otherwise set up the right-hand side for level i + 1, 

r i+1 G W m , ( r m, z i+i) = ( r Bi, EiZi+i) , Vz i+1 G Wm = U i+ i, (62) 

end. 

for i = L- 1,... ,1, 

Average the approximate corrections on substructure interfaces on level i, 

u B i = Ei (WAi + Ui+i) • (63) 

Compute the interior post- correction on level i, 

vu eUji-.a (vn, z u ) = a (u B i,zn) , Vz/* G Un. (64) 

Apply the combined corrections, 

Ui = un + u B i ~ vn. (65) 

end. 

We can now show that the Multilevel BDDC can be cast as the Multispace BDDC on energy 
orthogonal spaces, using the hierarchy of spaces from Figure [TJ 



Lemma 18 The Multilevel BDDC preconditioner in Algorithm\ll\is the abstract Multispace BDDC 
preconditioner from Algorithm^ with M = 2L — 1, and the spaces and operators 

X = U U Vi = U n , V 2 = (I - Pi)W A1 , V 3 = U I2 , 
V± = (I-P 2 )W A 2, V 5 = U I3 , ... (66) 
V2L-4 = (I - Pl-2)W A l-2, V 2L . 3 = t/ /L _i, 

V2L-2 = (I- P L -l)W AL -l, V 2L -1 = WuL-1, 

Qi=I, Q2 = Qs = (I-Pi)E 1 , 

Q 4 = Q 5 = (I-P 1 )Ei(I-P 2 )E 2 , ... (67) 

Q2L-4 = Q2L-3 = (/ - Pi) E 1 ■ ■ ■ (I - P L _ 2 ) E L _ 2 , 
Q2L-2 = Q2L-X = (I~Pl)E 1 ■■■ (I- P L _i) El-1, 

and the assumptions of Corollary [|] are satisfied. 
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Proof. Let r\ G U[. Define the vectors Vk, k = 1, . . . : 2L — 1 by ([3"|) with the spaces and operators 
given by (|66]) - ([67|) . and let un, r B i, WAi, wm, fi+i, ub%, vn, and Uj be the quantities in Algorithm H71 
defined by (|58 |> -(|B5 j> . 

First, with V\ = Un, the definition of v\ in (j3J) is ([58} with z = 1 and un = «i. We show that in 
general, for level i = 1, . . . , L — 1, and space k = 2i — 1, we get ([3]) with = so that Vk = un 
and in particular t>2L-3 = MiX-i- So, let z/j S Un, i = 2, . . . , L — 1, be arbitrary. From ([581) using 
d62D and 

a(uji,zn) = (ri,zu) = (rBi-i,Ei-izu) = (68) 
= (n-i, Ei-izu) - a (un-1, Ei-iZii) . 

Since from (|58|) using the fact that P{-\Ei_\zn & Un-i it follows that 

(ri_i, Pi-iEi^za) - a (uji-i, Pi-iEi-izn) = 0, 

we get from ([68]) . 

a(un, zji) = (n-i, (J - P_i) E^izn) - a (u/j-i, (J - P i _i).Ei_izj i ) , 

and because a (un-i, (I — Pi-i)Ei-\Zn) = by orthogonality, we get 

a(uu, zii) = (n-i, (I - Pi-i) Ei-xzji) ■ 

Repeating this process recursively using (|68p . we finally get 

a{un, zh) = (rj_i, (I - Pj_i) Ei^xZn) = ... 

= {n, (I — P\) Ei ■ ■ ■ (I — P^i) E^xzn) . 

Next, consider WAi € W A i defined by ([60|) , We show that for i = 1, . . . , L — 1, and k = 2i, we 
get © with Vk = WAi, so that v& = WAi and in particular V2L-2 = wal-i- So, let ZAi € WAi be 
arbitrary. From ([oTJj) using ([5U]) . 

a (io Ai , z A i) = (rBi,EiZ Ai ) = (r^EiZAi) ~ a (un, EiZ A i) ■ (69) 

From the definition of un by (|58p and since PiEiZAi € C/ji it follows that 

(ri,PiEiZAi) - a (un, PiEiZAi) = 0, 

so (j69[) gives 

a (ifAi, ZAi) = {n, (I - Pi) EiZAi) - a (un, (I - Pi) EiZAi) ■ 
Next, because a (un, (I — Pj) E^ZAi) = by orthogonality, and using ([62]) . 

a (wAi, ZAi) = {n, (I - Pi) EiZAi) = {rBi-i, E^ (I - Pi) EiZAi) - 
Repeating this process recursively, we finally get 

a(wAi,ZAi) = {ri,(I - Pi) EiZAi) = ■ ■■ 

= (n, (I -Pi) Ei ■■■ (I -Pi) EiZAi). 

To verify ([3]), it remains to show that PiWAi = 0; then WAi € (I — Pi)Waj = Vk- Since Pj is an 
a-orthogonal projection, it holds that 

a (PiWAi, PiWAi) = a(wAi,PiWAi) = (rBi, EiPiWAi) = 0, 
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where we have used EiUji C Un following the assumption ([57]) and the equality 

(rBh z n) = i r i, z n) ~ a ( u n, z n) = 

for any zn € Un, which follows from ([58]) and ([59]) , 

In exactly the same way, we get that if wul-i G Wul-i is defined by ([6T]) . then V2L-1 = wul-i 
satisfies ([3]) with k = 2L — 1. 

Finally, from ([63]) - ([65]) for any i = L — 2, . . . , 1, we get 

= un + (I - Pi) Ei (wAi + u i+ i) 

= un + (I - Pi) Ei [wAi + un+i + (I - Pi+i) E i+ i (wAi+i + u i+2 )] 
= un+ 

+ (I - Pi) Ei [w Ai + ... + (/- P L -i) El-\ {wal-\ + UnL-l)} , 
and, in particular for u%, 
ui = un + 

+ (J - Pi) Pi {WAX +... + (/- Pl-i) (WAL^I + UnL-l)] 

= Q\V\ + Q2V2 + • • • + Q2L-2V2L-2 + Q2L-1W2L-1- 

It remains to verify the assumptions of Corollary [4j 

The spaces Wm and Wa%, for all i = 1, . . . , L — 1, are a-orthogonal by ([33]) and from ([331) . 

(/-Pi) WAi C WAi, 

thus (I — Pi) WAi is a-orthogonal to Wni- Since Whi = Ci+i consists of discrete harmonic functions 
on level i from ([36]) . and Un+i C J7»+i, it follows by induction that the spaces given by ([66]) . 
are a-orthogonal. 

We now show that the operators Qk defined by ([67]) are projections. From our definitions, coarse 
degrees of freedom on substructuring level i (from which we construct the level i+1 problem) depend 
only on the values of degrees of freedom on the interface Tj and Tj C Ti for j >i. Then, 

(7 - P,-)£,-(Z - Pi)Ei{I - PfiEj = {I- Pi)Ei(I - PfiEj. (70) 

Using ([TO]) and since (I — Pi) Pi is a projection by ([3T]) . we get 

[(/ - Pi) Pi • • • (/ - P) P,] 2 = (/ - Pi) Pi (/ - Pi) Ex ■ ■ ■ (I - Pi) Pi 

= (I-P 1 )E 1 ---(I-P i )E i , 

so the operators Qk from ([67]) are projections. 

It remains to prove the decomposition of unity (jl4h . Let Ui € £/j, such that 

u'i = un + WAi + Uj+i, (71) 

u/j e 17/i, ifAi G (I - Pi) WAi, Ui+i G f/i+i (72) 

and 

«i = uji + (J - P) EiWAi + (I ~ Pi) EiU i+ x- (73) 
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From d?!]), WAi + G t^i since Uj G E/j and u/j 6 f/jj C Ui. Then i% [waj + u i+ i] = WAi + iti+i 
by ((57]), so 

fi = U/j + (I - Pi)Ei [WAi + = uu + Pi) [w A i + u i+1 ] = 

= Uu + WAi + = Uu + WAi + «t+l = *4, 

because wa; an d ^i+i are discrete harmonic on level i. The fact that Uj + i in (|TTj) and (|73"j) are the 
same on arbitrary level i can be proved in exactly the same way using induction and putting Uj+i 
in dnj as 

■Ui+i = nji + i + . . . + wal-1 + wni-i, 

G f/ji+i, wal-1 G (I - Pl-i) Wal-i, wul-i G Wul-i, 

and in ([73]) as 

u i+ i = un+i + ... + (/- Pj+i) E i+ i ■••(/- Pl_i) (-^al-i + wnx-i) • 

■ 

The following bound follows from writing of the Multilevel BDDC as Multispace BDDC in 
Lemma [181 and the estimate for Multispace BDDC in Corollary HI 

Lemma 19 If for some u > 1, 

\\(I-Pi)E lW Ai\\ 2 a <u;\\wAif a Vwai G (I-Pi)Wau 
||(J - Px)E x u I 2\^ a < U) ||w/2||o Vu/2 G U I2 , 

(74) 

||(7 - P X )E 1 •••(/- P L -i)E L -xw U l-i\\1 < oj ||x«iiL-i|la V^nL-i G Wul-i, 
then the Multilevel BDDC preconditioner (Algorithm\17\) satisfies k < u. 



Proof. Choose the spaces and operators as in (|66j) - (|67j) so that un = v% G V\ = Uu, 

WAi = V2 G V 2 = (I — Pi) Wai, Wul-1 = V2L-1 G V 2 L-i = Wul-1- The bound now follows 
from Corollary [H ■ 

Lemma 20 If for some Ui > 1, 

\\{I-Pi)EiWi\\l<uJi\\wi\\l, VwiGWi, i = l,...,L-l, (75) 

then the Multilevel BDDC preconditioner (Algorithm\17\) satisfies k < nfci 1 ^- 



Proof. Note from Lemma [191 that (I - Pi) Wai C Wai C Wi, U I2 C W m C Wi, and generally 
(I - Pi) W Al C Wa 1 C Wi, U Ii+ i c W m C Wi. ■ 

7 Condition Number Bound for the Model Problem 

Let 

l u; L(f2?) ^ e ^ ne energy norm of a function w G Wni, * — 1, ■ ■ ■ , L — 1, restricted to subdomain 
fif, s = 1, . . . Ni, i.e., |w|a(nf) = Jqs^wVw dx, and let ||n>|| a be the norm obtained by piecewise 
integration over each fif . To apply Lemma l20l to the model problem presented in Section [5j we need 
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to generalize the estimate from Theorem [TBI to coarse levels. To this end, let : Wm — > C^i+i be 
an interpolation from the level i coarse degrees of freedom (i.e., level % + 1 degrees of freedom) to 
functions in another space U^x and assume that, for all i = 1, . . . ,L — 1, and s = 1, . . . , Ni, the 
interpolation satisfies for all w G Wrij and for all the equivalence 

ci,i |I i+ iu;^ (nf+i) < |/i^|^ ( n f+1 ) < <H,2 l J m^la(nf +1 ) > ( 76 ) 
which implies by Lemma [22] also the equivalence 

Q,l |^+lw|^l/ 2(9 ^ +i) < |-?iHtfV2(Sfif +1 ) ^ C i,2 l^+l^l/fl/a^Q^J , (77) 

with Ci^/cj.i < const bounded independently of ilo, ■ ■ • , Hi+%. 

Remark 21 Since I\ = I, the two norms are the same on Wno = U\ = U\. 

For the three-level BDDC in two dimensions, the result of Tu [201 Lemma 4.2], which is based on 
the lower bound estimates by Brenner and Sung [2], can be written in our settings for all w € Wm 
and for all 0| as 

C1.1 I^Ha(ng) < kla(Q|) ^ C l,2 I^Ha^f) , (78) 

where -Z2 is a piecewise (bi) linear interpolation given by values at corners of level 1 substructures, 
and ci 5 2/ci 5 i < const independently of H/h. 

For the three-level BDDC in three dimensions, the result of Tu |19[ Lemma 4.5], which is based 
on the lower bound estimates by Brenner and He [TJ, can be written in our settings for all w £ Wm 
and for all 0| as 

ci,i |£H#i/2(an|) ^ \ w \m/2(dn°) ^ c i,2 Ihwfm^idQ^ > ( 79 ) 

where I2 is an interpolation from the coarse degrees of freedom given by the averages over 
substructure edges, and cx^/cia < const independently of H/h. 

We note that the level 2 substructures are called subregions in [2Ql [19] . Since I\ = I, with i = 1 
the equivalence ([75]) corresponds to ([76]) . and ([79]) to ([77]) . 

The next Lemma establishes the equivalence of seminorms on a factor space from the equivalence 
of norms on the original space. Let V C U be finite dimensional spaces and ||-|L a norm on U and 
define 

\ u \u/v.A = m jn\\u-v\\ A . (80) 

We will be using (|80p for the norm on the space of discrete harmonic functions (/ — P%)Wi with 
V as the space of interior functions Un, and also with V as the space Wai- In particular, since 
Wn; C (I - Pi)Wi, we have 

w € Wn;, HHL = mi ii ll^-^Alla (81) 

Lemma 22 Xei IHIb be norms on U, and 

ci IMIs ^ ^ c 2 IMIs > Vu £ U. (82) 



T/ien /or any subspace V C U , 



It 2 ^ I |2 ^ I t2 

c l l 7 %/V,.B ^ \ U \U/V,A - c 2 l n l(7/V,B ' 



resp., 

ci min ||u — v\\* B < min |ju — < C2 min ||u — v\\* B . (83) 
v&V v£V vGV 
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Proof. From the definition (|80p of the norm on a factor space, we get 

Hu/V,A = \\ u ~ v \\a = \\ u ~va\\ a . 
for some va- Let vb be defined similarly. Then 



\u\ulVA = min ll n ~~ v \\a = W u ~ v a\\a ^ \\ u ~ v b\\a ^ 



< c 2 \\u - vb\\ 2 b = C2min ||« - v\\ B = c 2 Mjy/y,B 



which is the right hand side inequality in (|83|) . The left hand side inequality follows by switching 
the notation for ||-||^ and \\-\\b- ■ 



Lemma 23 For all i = 0, . . . , L — 1, and s = 1, . . . , N t , 

|2 ■ . 1 2 



Ci,i |Ii+H a(n|+i) < M a (nf +1 ) < c fe,2 |Ji+i^L(of +1 ) » Vw e ^m, VO| +1 , (84) 
wrf/i Cj^/cj,! < Cj> independently of Hq,. . . , Hi + \. 



Proof. The proof follows by induction. For i = 0, ([84]) holds because I\ = I. Suppose that (|84p 
holds for some i < L — 2 and let u; G Wni+i- From the definition of Wni+i by energy minimization, 

\ w \ a (n° i+1 ) = "il 11 \w-w A \ a{n s ). (85) 

From (|85p . the induction assumption, and Lemma [22] eq. (|83p . it follows that 

(HA min [Ji + i«; - Ji+iiOAlofn* ) ( 86 ) 
< min |w-u> A |a(c- ) < c i,2 min \I l+1 w - I i+1 w A \ 2 ain s ) 

From the assumption (|76j) . applied to the functions of the form Ii+\w on 2 , 

Cj \Ii + 2w\ 2 a(af j < min |J i+1 u; - J i+1 u; A £ (n? } < c 2 \I l+2 w\ 2 a(Q3 j (87) 

with c 2 /ci, bounded independently of H ,...,H i+1 . Then ([55]) . (f86|) and ([57]) imply ([51]) with 
Q = Cj_ic 2 /ci. ■ 

Next, we generalize the estimate from Theorem [15] to coarse levels. 

Lemma 24 For sub structuring levels i = 1, . . . , L — 1, 

\\(I-P i )E i w i \\l<C i (l + ]og^^\ \\wi\\l, y Wi eUi. 
Proof. From (|84p . summation over substructures on level i gives 

c *,l ll-^HIa - HHIa - C «>2 ll-^HIa ) G f/j. (89) 

Next, in our context, using the definition of Pj and (|83|) . we get 

- = min WhEiWi - IiUn\\l , 

uueUu 

20 



so from (J52J) for some Ci and all i = 1, . . . , L — 1, 

WhEiWi - IiUnWl < d I 1 4- Iok „ * 
Similarly, from ([8"9"j) and (|9"0"j) it follows that 



min WhEiWi - IiUn\\ 2 a < d ( 1 + log ) 11^^11^ , Vu>j G f/j. (90) 



which is (I88H with 



|(I - P i )E i w i \\ a = min \\E i w i -u I i\\ a 

<c i<2 min | \IiEiWi - huu ||^ 

< Ci^Ci ^1 + log "^-^ ll^illa 



L>i — , 

CiA 



and c^/c^i from Lemma 

Theorem 25 The Multilevel BDDC for the model problem and corner coarse function in 2D and 
edge coarse functions in 3D satisfies the condition number estimate 



^aiYc^i+iog 



m ^ 2 



Proof. The proof follows from Lemmas [20] and [231 with uji = Ci ^1 + log J Ii ^ 

Remark 26 For L = 3 in two and three dimensions we recover the estimates by Tu \2C\ [7 
respectively. 

Remark 27 While for standard (two-level) BDDC it is immediate that increasing the coarse space 
and thus decreasing the space W cannot increase the condition number bound, this is an open 
problem for the multilevel method. In fact, the 3D numerical results in the next section suggest that 
this may not be the case. 

Corollary 28 In the case of uniform coarsening, i.e. with Hi/Hi_i = H/h and the same geometry 
of decomposition on all levels i = 1, . . . L — 1, we get 

k < C L - 1 (1 + log H/h) 2{L - 1] . (91) 



8 Numerical Examples 

Numerical examples are presented in this section for the Poisson equation in two and three 
dimensions. The problem domain in 2D (3D) is the unit square (cube), and standard bilinear 
(trilinear) finite elements are used for the discretization. The substructures at each level are squares 
or cubes, and periodic essential boundary conditions are applied to the boundary of the domain. 
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Figure 2: Plot of data in Tables [T] and [5] for Hi/H^i = 3 at all levels. 
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Figure 3: Plot of data in Tables [1] and [2] for uniform coarsening 
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Table 1: Two dimensional (2D) results. The letters C and E designate the use of corners and edges, 
respectively, in the coarse space. 



L 


C 


C+E 


n 






iter 


cond 


iter 


cond 










Hi/Hi- 


i = 3 at all levels. 




2 


8 


1.92 


5 


1.08 


144 


80 


3 


13 


3.10 


7 


1.34 


1296 


720 


4 


17 


5.31 


9 


1.60 


11,664 


6480 


5 


23 


9.22 


10 


1.85 


104,976 


58,320 


6 


31 


16.07 


11 


2.12 


944,748 


524,880 


7 


42 


28.02 


13 


2.45 


8,503,056 


4,723,920 






Hi/Hi- 


i = 4 at all levels. 




2 


9 


2.20 


6 


1.14 


256 


112 


3 


15 


4.02 


8 


1.51 


4096 


1792 


4 


21 


7.77 


10 


1.88 


65,536 


28,672 


5 


30 


15.2 


12 


2.24 


1,048,576 


/ICO 7CO 

458,752 





42 


29.7 


13 


2.64 


16,777,216 


i ,3411,032 






Hi/Hi- 


i = 8 at all levels. 




2 


10 


2.99 


7 


1.33 


1024 


240 


Q 
O 


19 


7.30 


11 


2.03 


65,536 


lo,oOU 


4 


31 


18.6 


13 


2.72 


4,194,304 


983,040 


5 


50 


47.38 


15 


3.40 


268,435,456 


62,914,560 






Hi/ Hi- ] 


= 12 at all levels. 




2 


11 


3.52 


8 


1.46 


2304 


368 


3 


21 


10.12 


12 


2.39 


331,776 


52,992 


4 


39 


29.93 


15 


3.32 


47,775,744 


7,630,848 






Hi/ Hi- ] 


= 16 at all levels. 




2 


11 


3.94 


8 


1.56 


4096 


496 


3 


23 


12.62 


13 


2.67 


1,048,576 


126,976 


4 


43 


41.43 


16 


3.78 


268,435,456 


32,505,856 



This choice of boundary conditions allows us to solve very large problems on a single processor 
since all substructure matrices are identical for a given level. 

The preconditioned conjugate gradient algorithm is used to solve the associated linear systems 
to a relative residual tolerance of 10 -8 for random right-hand-sides with zero mean value. The zero 
mean condition is required since, for periodic boundary conditions, the null space of the coefficient 
matrix is the unit vector. The coarse problem always has 4 2 (4 3 ) subdomains at the coarsest level 
for 2D (3D) problems. 

The number of levels (L), the number of iterations (iter), and condition number estimates (cond) 
obtained from the conjugate gradient iterations are reported in Tables [1] and [2j The letters C, E, 
and F designate the use of corners, edges, or faces in the coarse space. For example, C+E means 
that both corners and edges are used in the coarse space. For 2D and 3D problems, the theory is 
applicable to coarse spaces C and E, respectively. Also shown in the tables are the total number 
of unknowns (n) and the number of unknowns (nr) on subdomain boundaries at the finest level. 

The results in Tables Q] and [2] are displayed in Figure [2] for a fixed value of Hi/Hi_i = 3. In 
two dimensions we observe very different behavior depending on the particular form of the coarse 
space. If only corners are used in 2D, then there is very rapid growth of the condition number with 
increasing numbers of levels as predicted by the theory. In contrast, if both corners and edges are 
used in the 2D coarse space, then the condition number appears to vary linearly with L for the the 
number of levels considered. Our explanation is that a bound similar to Theorem [25] still applies to 
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Table 2: Three dimensional (3D) results. The letters C, E, and F designate the use of corners, 
edges, and faces, respectively, in the coarse space. 



T 

L 


E 


C+E 


C+E+F 


n 


nr 




iter 


cond 


iter 


cond 


iter 


cond 












i — 3 at all levels. 






2 


10 




1.85 


8 


1.47 


5 


1.08 


1728 


1216 


Q 
O 


14 




3.02 


12 


2.34 


8 


1.50 


46,656 


32,832 


4 


18 




4.74 


18 


5.21 


11 


o on 


1,259,712 


886,464 





23 




7.40 


26 


14.0 


16 


3.98 


34,012,224 


23,934,528 








l = 4 at all levels. 






2 


10 




1.94 


9 


1.66 


6 


1.16 


4096 


2368 


3 


15 




3.51 


14 


3.24 


10 


1.93 


262,144 


151,552 


4 


20 




6.09 


22 


9.95 


14 


3.05 


16,777,216 


9,699,328 






Hi/ Hi- 


i = 8 at all levels. 






2 


12 




2.37 


11 


2.24 


8 


1.50 


32,768 


10,816 


3 


19 




5.48 


20 


7.59 


14 


3.32 


16,777,216 


5,537,792 






Hi/ Hi-i 


= 10 at all levels. 






2 


12 




2.56 


12 


2.47 


9 


1.69 


64,000 


17,344 


3 


20 




6.39 


22 


10.1 


16 


3.85 


64,000,000 


17,344,000 



the favorable 2D case, though possibly with (much) smaller constants, so the exponential growth 
of the condition number is no longer apparent. The results in Tables [1] and [2] are also displayed 
in Figure [3] for fixed numbers of levels. The observed growth of condition numbers for the case of 
uniform coarsening is consistent with the estimate in (|9ip . 

Similar trends are present in 3D, but the beneficial effects of using more enriched coarse spaces 
are much less pronounced. Interestingly, when comparing the use of edges only (E) with corners and 
edges (C+E) in the coarse space, the latter does not always lead to smaller numbers of iterations or 
condition numbers for more than two levels. The fully enriched coarse space (C+E+F), however, 
does give the best results in terms of iterations and condition numbers. It should be noted that the 
present 3D theory in Theorem 1251 covers only the use of the edges only, and the present theory does 
not guarantee that the condition number (or even its bound) decrease with increasing the coarse 
space (Remark |27|) . 

In summary, the numerical examples suggest that better performance, especially in 2D, can be 
obtained when using a fully enriched coarse space. Doing so does not incur a large computational 
expense since there is never the need to solve a large coarse problem exactly with the multilevel 
approach. Finally, we note that a large number of levels is not required to solve very large problems. 
For example, the number of unknowns in 3D for a 4-level method with a coarsening ratio of 
H./Hi^ = 10 at all levels is (10 4 ) 3 = 10 12 . 
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